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Abstract 

We have used Asymptotic Iteration Method (AIM) for obtaining the eigenvalues 
of the Schrodinger’s equation for a deformed well problem representing trigonometric 
functions. By solving the problem, we have found that the Schrodinger’s equation 
for the considered potential has quasi-exact solutions. Additionally, we have also 
calculated the perturbation expansion of energy eigenvalues and found very simple 
analytical expression of the energy. Finally, we have considered more general cases 
and obtained energy eigenvalues for arbitrary potential parameters. 

Key Words: Asymptotic iteration method, quasi-exact solution, perturbation method, 
approximate solution 


1 Introduction 


The solution method of the Schrodinger’s equation, which is the main scope of quantum 
mechanics, depends upon the configuration of the system held, i.e., behavior of the po¬ 
tential. One of the advantages of the asymptotic iteration method (AIM), introduced in 
[l, 131, is that it can be used to solve the Schrodinger’s equation exactly |3. ffl R [fill 17f]] 
or approximetly for numerous potentials 


10|, 111, UJZ]- 


AIM can also be used in the framework of perturbation theory. In common perturbation 
theory, eigenfunctions of unperturbed hamiltonian (i.e. exact solvable part) is required 
to construct the perturbation expansion for energy eigenvalues and eigenfunctions. Fur¬ 
thermore, it is big challenging even for third or fourth-order of approximation in some 
cases. The advantage of the perturbation method using AIM is that the coefficients in 
perturbation expansion can be calculated directly without any need to eigenfunctions of 


unperturbed hamiltonian [5|, [9K1J 
We use the AIM to find out both the quasi-exact and approximate (numeric) solutions 
of the Schrodinger’s equation for a deformed well problem representing trigonometric 
functions in the present work. The AIM is summarized in section II. We find out the 
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eigenvalues of Schrodinger’s equation approximately and quasi-exactly in section III for 
a potential defined by 

V(x) = A Cos(x) + 5 2 Sin 2 (x) + x€( 0 , 7 r) (1) 

Sin (x) 

where A and B are coefficients and 7 ranges over (0, 00 ). In section IV, we apply 
perturbation method for this kind of potential , just after, summarizing the perturbation 
theory in the framework of AIM. 

2 The Asymptotic Iteration Method (AIM) 

Although the asymptotic iteration method has been introduced comprehensively in [1], 
we outline this method in this section according to organization of the paper. The AIM 
can be used to solve second-order homogeneous linear differential equations of the form 

y"{x) = A 0 (x)y'(x) + s 0 (x)y(x) ( 2 ) 


where Ao and so have continuous derivatives in the defined interval of the x independent 
variable. According to asymptotic aspect of the method, if we have 


Sn 

Sn-1 


A 


n 


A 


71—1 


= a 


(3) 


for sufficient large n > 0. The general solution of Eq{2] is obtained as 


y(x) = exp — / a(t)dt 


C 2 + C 1 / exp / (A 0 (t) + 2a{t)) dr dt 


(4) 


where 


X n — A 'n-\ + Sn-l + AoA ra -l (5) 

Sn — ^n —1 T SoAn—1 

Eigenvalues, E, of the Schrodinger’s equation are obtained using the following termina¬ 
tion condition by means of Eql3] 

S n (x) = s n (x)X n -i(x) - X n (x)s n -i(x) = 0 ( 6 ) 


where the eigenvalues of n energy levels are obtained by Eql 6 ]with n iterations, if the 
problem is analytically solvable (exact solution). In this case, eigenvalues and eigen¬ 
functions can be determined in explicit algebraic form in AIM. Furthermore, there are 
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limited number of suitable potentials for this case. 

In exact solvable problem, condition 5 n = 0 is satisfied at every x point in the defined 
interval since 5 n is independent from x. The (unknown) E values are determined from 
the n roots of condition 5 n = 0. In case of approximation, 5 n depends on both x and 
E. So, we have to determine a, x = xq initial value for solving 5 n = 0 with respect to E. 
As this initial xq value may be found from the minimum value of the potential may also 
bounded domain of the eigenvalue problem a Q. 


3 AIM for the Problem 

We have considered a deformed well problem representing trigonometric functions as 
below 


V(x) = A Cos(x ) + i? 2 Sin 2 (x) + 


2 0 ; „ 2 ,^ , 7(7 + 1) 

2 ... ’ 


Sin 2 x 


x € ( 0 , 7 r) 


where A, B and 7 are the potential parameters. 7 is defined in the interval [0, 00 ). If 
we substitute V (x) in the Schrodinger’s equation, we have 


f d 2 ^ 
\i^ + E ~ 


2 o :_ 2 /_n , 7(7 + 1) 


A Cos(x) + B 2 Sin 2 (x) + 


Sin 2 (x) 


'k(x) = 0 


(7) 


where unnormalized eigenfunction ’I'(x) € L 2 (0,ir). We may assume that T(x) takes 
the form 


+ (x) = ShP +1 (x) e BCo ^ x) f{x) 


( 8 ) 


with this choice, T(x) vanishes at x = 0 and x = ir. Once we substitute above equation 
in EqU] 


/"(*)+2 


( 7 + 1 ) i^~ BSin(x) 


f'(x)+ E - (7 + l ) 2 - {A + B (27 + 3)} Cos(x) f(x) = 0 

(9) 


is obtained. If we change the variable such as Cos(x) = y in Eqj9] which yields 

f "( y ) = 2 (1^5 - B ) /'to) + (yt^) /(») ( 10 ) 

where uj = E — (7 + I) 2 , £ = A + 2aB and 2 a = 27 + 3 . The last equation is in form 
of Eq{2] for which AIM can be used. We obtain the (unkown) E energy eigenvalues by 
using AIM. 

3.1 Quasi-Exact Solutions of the Problem 

We have obtained Eq(T0] for using the AIM for determining the E eigenvalues of EqJTl 
Our calculations show that the problem has quasi-exact solution under some special 
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conditions on the potential parameters. AIM gives us that if £ = —2nB with n = 
0,1, 2,3,... the energy eigenvalues can be obtained exactly. We will give first two results 
analytically because of the complecity of the calculations, but we will give more numerical 
results in Table]!] Using £ = —2 nB and £ = A + 2aB together, we get the following 
formula 


A =-2B (^n + 'y+^j (11) 

For n = 0, we get oj = 0 and Eqq = (7 + l) 2 . For n = 1, we get lo = a- F V + 4.B 2 and 
the energies can be obtained as follows 

= (y+£f+ir\/(j +£) 2+4B2 < 12 > 

Eu = (^ + 5) 2 + i + \/(^ + 5) 2 + 4B2 (13) 

Numerically calculated eigenvalues for B = 1 and a = 4.5 can be found in Tabled 


Table 1. Numerically calculated eigenvalues for B = 1 and a = 4.5 for n = 0 to 5 by 
using £ = —2 nB. k is sub-level for each n level and k= 0, 1, 2, ..., n. 


n 

k 

A 

^nk 

Enk 

0 

0 

-11 

0 

16.0000 

1 

0 

-13 

-0.424429 

15.5756 


1 

-13 

9.42443 

25.4244 


0 

-15 

-0.919071 

15.0809 

2 

1 

-15 

9.27711 

25.2771 


2 

-15 

20.642 

36.642 


0 

-17 

-1.47815 

14.5218 

3 

1 

-17 

9.09552 

25.0955 


2 

-17 

20.6205 

36.6205 


3 

-17 

33.7621 

49.7621 


0 

-19 

-2.09614 

13.9039 


1 

-19 

8.87768 

24.8777 

4 

2 

-19 

20.5892 

36.5892 


3 

-19 

33.7953 

49.7953 


4 

-19 

48.834 

64.834 


0 

-21 

-2.76786 

13.2321 


1 

-21 

8.62229 

24.6223 


2 

-21 

20.5459 

36.5459 

5 

3 

-21 

33.8293 

49.8293 


4 

-21 

48.8905 

64.8905 


5 

-21 

65.8799 

81.8799 


3.2 Approximately Determination of the Eigenvalues 

In this case, we have determined the eigenvalues for arbitrary A, B and 7 values. Cal¬ 
culated eigenvalues (E c i ire ,ct) of n = 0 to 5 are given in Tabled in which A, B and 7 (in 
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(A; B ; 7 ) form) are (0.167; 0.0019; 0.0019), (2; 0.5; 1) and (4; 1; 1.97). In the calculation 
we take yo = k- 


Table 2. Approximately calculated eigenvalues Edirect for (0.167; 0.0019; 0.0019), (2; 
0.5; 1) and (4; 1; 1.97) in (A; B ; 7 ) form. The value written as subscript for each 
eigenvalue represents the iteration number. 



A = 0.167, B = 7 = 0.0019 

A = 2 , B = 0.5, 7 = 1 

A = 4, B = 1 , 7 = 1.97 

n 

LO 

direct 

LO 

direct 

LO 

direct 

0 

-0.00231 

1.00149(3) 

0.07399 

4.07399 a o) 

0.57328 

9.39418(13) 

1 

3.00473 

4.00853(4) 

5.17104 

9.17104 a2 ) 

7.618 

16.4389(13) 

2 

8.008 

9.0118( 5 ) 

12.1609 

16.1609(11) 

16.5264 

25.3473(i 4 ) 

3 

15.0116 

16.0154(g) 

21.1506 

25.1506(13) 

27.4331 

36.254(16) 

4 

24.0153 

25.0191 m 

32.1438 

36.1438(14) 

40.348 

49.1689(18) 

5 

35.0191 

36.0229(g) 

45.1392 

49.1392 ( i 5 ) 

55.2698 

64.0907(i 9 ) 


As is seen from Tabled energy eigenvalues act as ~ (n + 7 + l) 2 . 

4 Perturbation Theory within the frame of AIM 

Although the usage of perturbation method within AIM is introduced in [ 9 ] comprehen¬ 
sively, we give a summary about the application method in this section. Assume that 
the potential of a system is written in a form of 

V(x) = V 0 (x) + nV p (x) (14) 


where Vq(x) is solvable (or unperturbed hamiltonian) potential. V p (x ) and /r are potential 
of the perturbed hamiltonian and perturbation expansion parameter, respectively. So, 
the Schrodinger’s equation reads 

+ V °( X ) + V V P( X )) *( X ) = E *( X ) ( 15 ) 

where E n eigenvalues are written as the expansion of j th-order correction Eft as follows 

OO 

En = 4 0) + ^ + y 2 E^ + ... = jr (16) 

j =0 

The form of EqJ2]can be yielded by a substitution of (in general form) ip{x) = ip 0 (x)f(x) 
in Eall51 So, one can obtain the following equation for f(x) 

f"(x) = A 0 (x, /i, E)f'(x) + s 0 (x, n, E)f(x) (17) 


The termination condition, EqEl is written as 

5 n (x, fi, E) = s n (x, n, E)X n -i(x, n, E) - X n (x, n, E)s n - i(x, y, E) = 0 (18) 
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Once 5 n (x, fi, E) is expanded about n = 0, we obtain 


r /_ .. tti\ r n m ^ d8 n {x,H,E) 

o n \% i l-ii E) — o n (x, 0 , .£/)+ ^ 


/i 2 d 2 5 n (x,iJ,,E) 


/i=0 


2 ! 


0/r 2 


/i=0 


+... = J2» k 6W(x,E) = 0 

(19) 


fc =0 


where 8n\x,E ) 


1 d k Sn(x,/l,E) 
k\ 


fi=0 


The condition of 


<f)(x,£) = 0 


( 20 ) 


can be obtained by means of Eq(19]for each k. According to the perturbation method in 
the framework of AIM, solving the equation 5 n (x, 0 ,E) = 0 with respect to (unknown) E 
gives En ^ (eigenvalues of unperturbed hamiltonian), equation 5n\x,E) = 0 gives E^ 
(first-order correction to E n ), 5n (x,E) gives En ; (second-order correction to E n ) and 
so on. The eigenfunctions can also be found with similar manner. This is an attractive 
feature of the AIM usage in the perturbation theory for obtaining the eigenfunctions 
f n {x) given as follows 


f n (x) = exp 


a n (t, / i)dt 


( 21 ) 


where a n (t, fi) = s n (t, n)/\ n (t, fi). a n (t , /r) is expanded about /r = 0 in a similar manner 
of obtaining the eigenvalues. So, 


a n (t,n) = ^(j k a^\t) 
k=0 


( 22 ) 


(k) 

where a\> x) 
follows 


1 & ° i n(x,ii) 

H dp 


/x=o 


Thus, perturbation expansion of the f n t x \ is written as 


f n (x) = exp 




II fnH*) 

k=0 


where k th-order correction 


fn\x) to fn(x) is 



(23) 


( 24 ) 
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4.1 Perturbation Theory for the Potential Function in the frame of 
AIM 

In section III, we have obtained the eigenvalue equation in the appropriate form for AIM 
usage as follows 

/"<#) = 2 - B ) /'(») + (f 'Ajf) f(v) ( 25 ) 

where u ; = E — (7 + l) 2 , £ = A + 2aE and 2a = 27 + 3. A and B parameters of 
the potential V(x), given by Eq(H was linked to eachother via EqQT] Suppose that 
A = 2 aB, in this case £ = 2 (a + a)B and the differential equation can be written as 
below 


^) = 2 (r?- 8 )^ + ( 2(l+ .-M ^ <*> 

Now, B can be taken as perturbation expansion parameter. We can expand oj as below 

u n = uW + BuV + B 2 uW + ... (27) 


In this expansion the general form of zeroth-order correction oj 


(o) 


is obtained via, 


^ 0 )(x,0,u4 0) )=0 


(28) 


as mentioned previously. The obtained few values of w® are 0 , 2a, 2 + 4a, 6 + 6 a, 
12 + 8 a, 20 + 10 a, respectively. One can easily write general formula for ojn as below 

uj$ = n (n — 1) + 2 na (29) 


The first-order correction to ui n is generalized in the similar manner, and it is found 
that ujn J = 0 for each n level via the equation <5^ (a:, 0, ) = 0. The second-order 

correction to oj n is generalized as 


( 2 ) _ 0 f a 2 gi{n,a) - g 2 (n,cr)\ 

- ) 


(30) 


where 


gi (n, a) = 2a 2 — (2 n + 5) a — (n 2 — n — 3) 

g 2 (n, a) = 2a 1 + (6 n — 5) a 3 + {in 2 — 11 n + 3) a 2 + 4n (n — l) 2 a + n 2 ( n — 1)^31) 

53 (n, a) = (2a + 2n + 1) (2a + 2n — 3) (a + n) (a + n — 1) 
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If Eq|29] and Eqj30] are substituted in Eqj27] one can obtain 


u n = n (n — 1) + 2ncr 


A 2 gi{n,a) - AB 2 g 2 (n,a) 
253(n,cr) 


( 32 ) 


where we use a = substituted in Eq|30]as well. From u n = E n — (7 + l) 2 , we get En 
(perturbed energy) approximately as follows 


E^( 7 + n + l) 2 


A 2 gi(n,a) - AB 2 g 2 {n,a) 
2ff3(n,cr) 


(33) 


Numeric results are given in Tablc|3]in which En, given in Eg 1331 was used for n = 0 to 
5, and in Tablejd] for n = 0 to 10. The comparison of En with numerically calculated 
eigenvalues E^irect for arbitrary A, B and 7 is also given in Tablej3] and TablcJdJ The 
values written as subscript represent the iteration numbers in both tables. 


Table 3. Comparison of eigenvalues of perturbed hamiltonian En with numerically ones 
Edirect for A = 2, B = 2 and 7 = 4 for n = 0 to 5. 


n 

^ n 

direct 

E P n 

0 

3.6259 

28.6259(13) 

28.6364 

1 

14.1219 

39.1219( 14 ) 

39.1329 

2 

26.8245 

51.8245( 17 ) 

51.8308 

3 

41.6323 

66.6323(18) 

66.6353 

4 

58.5005 

83.5005(i 9 ) 

83.5015 

5 

77.406 

102.406(i 9 ) 

102.406 


Table 4. Comparison of the eigenvalues of perturbed hamiltonian En with numerically 
ones Edirect for A = 2, B = 1 and 7 = 0.0019 for n = 0 to 10. 


n 

Wn 

Edirect 

En 

0 

0.38637 

1.39017(i 2 ) 

1.42145 

1 

3.65921 

4.66301(13) 

4.64073 

2 

8.56603 

9.56983(i 6 ) 

9.56855 

3 

15.5445 

16.5483(i 5 ) 

16.547 

4 

24.5364 

25.5402(i 6 ) 

25.5392 

5 

35.5338 

36.5376(i 8 ) 

36.5368 

6 

48.5336 

49.5374(18) 

49.5369 

7 

63.5349 

64.5387(19) 

64.5383 

8 

80.537 

81.5408( 2 o) 

81.5404 

9 

99.5392 

100.543(20) 

100.543 

10 

120.542 

121.546( 2 i) 

121.546 


It can be seen from Tabic 13] that numerically calculated eigenvalues Edi rec t are consistent 
with the ones of perturbed hamiltonian En- Also, they are more consistent as the energy 
level n increases, and the values of Edirect converge more slowly with increasing n. All 
these comments apply to TableQ]as well. Futhermore, the eigenvalues act as ~ (n+ 7 +l ) 2 
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while oj n ~ n(n + 2) in Table|4j 

In TableEl values of E^ rect and E p x are inconsonant by selecting A = B = 5 and 7 = 0.5. 

Table 5. Comparison of the eigenvalues of perturbed hamiltonian E p with numerically 
calculated ones Edirect for A = B = 5 and 7 = 0.5 for n = 0 to 5. The values written as 
subscript represent the iteration numbers. 


n 

^ n 

E direct 

E p n 

0 

11.8935 

14.1435(3!) 

21.0000 

1 

19.3048 

21.55 48(34) 

20.8333 

2 

25.0743 

27.3243(3 6 ) 

25.7917 

3 

32.0111 

34.2611(40) 

33.375 

4 

41.5015 

43.7515(3 7 ) 

43.1667 

5 

53.2222 

55.4722(3 9 ) 

55.0476 


5 Conclusion 

We have used AIM for solving Schrodinger’s equation for a deformed well problem rep¬ 
resenting trigonometric function. In this study, we have applied AIM to quasi-exact 
solvable eigenvalue problem for obtaining the energy eigenvalues and seen that each en¬ 
ergy level n has a sub-level k that k = 0,1, 2,..., n. Method has also been used to obtain 
the eigenvalues approximately by using arbitrary potential parameters. According to re¬ 
sults, we can say that the energy eigenvalues behave as ~ (n + 7 + l ) 2 for larger quantum 
number n, as is seen in Eqj33j AIM can also be used in the framework of perturbation 
theory. In common perturbation theory, eigenfunctions of unperturbed hamiltonian is 
required to construct the perturbation expansion for energy eigenvalues and eigenfunc¬ 
tions. Also, determining the third or fourth-order correction to energy eigenvalues is 
usually quite big challenging. The advantage of the perturbation method using AIM is 
that the coefficients in perturbation expansion can be calculated directly without any 
need to eigenfunctions of unperturbed hamiltonian and any complication. 

Perturbation method has also been applied to potential function held in scope of the 
AIM, in this study. According to results, eigenvalues of perturbed hamiltonian are con¬ 
sistent with those of approximately calculated ones for which arbitrary parameters of the 
potential have been used. However, they are inconsonant for some set of the parameters. 
Although this inconsistency is true for n = 0 to 5, it tends to decrease as n increases. 
Consequently, it can be said that we can make the eigenvalues tend to specific behavior 
by selecting appropriate potential parameters. 
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